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ABSTRACT 

Context. The detection of trends or gradients in the transmission spectrum of extrasolar planets is possible with observations at very 
low spectral resolution. Transit measurements of sufficient accuracy using selected broad-band filters allow for an initial characteri¬ 
zation of the atmosphere of the planet. 

Aims. We want to investigate the atmosphere of the hot Jupiter HAT-P-12b for an increased absorption at the very blue wavelength 
regions caused by scattering. Eurthermore, we aim for a refinement of the transit parameters and the orbital ephemeris. 

Methods. We obtained time series photometry of 20 transit events and analyzed them homogeneously, along with eight light curves 
obtained from the literature. In total, the light curves span a range from 0.35 to 1.25 microns. During two observing seasons over four 
months each, we monitored the host star to constrain the potential influence of starspots on the derived transit parameters. 

Results. We rule out the presence of a Rayleigh slope extending over the entire optical wavelength range, a flat spectrum is favored 
for HAT-P-12b with respect to a cloud-free atmosphere model spectrum. A potential cause of such gray absorption is the presence of 
a cloud layer at the probed latitudes. Eurthermore, in this work we refine the transit parameters, the ephemeris and perform a TTV 
analysis in which we found no indication for an unseen companion. The host star showed a mild non-periodic variability of up to 1%. 
However, no stellar rotation period could be detected to high confidence. 

Key words. Techniques: photometric - Planets and satellites: atmospheres - Stars: individual: HAT-P-12 


1. Introduction 


Transit measurements at multiple wavelengths allow for measur¬ 
ing the so-called transmission spectrum of an extrasolar planet. It 
consists of the search for a wavelength dependence of the appar¬ 
ent radius of an exoplanet. At wavelengths of high atmospheric 
opacity, the atmosphere becomes optically thick at higher alti¬ 
tudes for the transmitted star light than at wavelengths of lower 
opacity. Such measurements are not limited to spectroscopic 
transit observation. In fact, spectrophotometry of transits can be 
done using multiple broad-band filters to create a “spectrum” 
that consists of a few data points over the optical and near in¬ 
frared (NIR) wavelength range. Similar to the color index of 
stars, which allows for a rough classification of stars, it might 
be possible to distinguish between several categories of plane¬ 


tary at mospheres by using only a handful well-se parated color 
filters (iNascimbeni et al.ll2013b[lBento et al.ll2014t) . 

The required accuracy when measuring the relative differ¬ 
ences in the apparent planetary radius can be achieved on small 
telescopes by using broad-band filters. Stromgren or even nar¬ 
rower filter systems in combination with a small telescope aper¬ 
ture would yield too high a photon noise in the light curves. 
However, the use of broad filters like the Johnson/Bessel or 
Sloan filter system comes with the drawback that signals of nar¬ 
row spectral features in the planetary atmosphere are strongly 
damped and smeared out. For example, detecting absorption 
from the sodium D line becomes more difficult. Instead, broad¬ 
band observations can be used to search for general trends in 
the spectrum, such as the rising absorption toward short wave¬ 
lengths caused by Rayleigh-scattering or the predicted redward 
slope between very blue colors and the central part of the optical 
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window caused by TiO (iFortnev et alJl201(ih . Despite the grow¬ 
ing number of observational publications in very recent years, 
the task remains challenging because of the weakness of the ex¬ 
pected signals. 

The atmospheric properties of several exoplanets have 
been constrained with gr ound- based broad- b and o bservations. 
iNascimbeni et al.l (l2013bh and iBiddle et al.l (12014l) show that 
the spectrum of the warm Neptune GJ3470b cannot be ex¬ 
plained with a gray, featureless atmosphere. Instead, a Rayleigh 
slope needs to be included in the modeling to gain a rea¬ 
sonable fit to the obser ved data. For t he at mosphere of 
the super-Earth GJ1214b, Ide Mooii et al.l (12012h rule d out a 
cloud-free solar composition, and INascimbeni et all (12015h 
provide d evidence agains t Rayle igh scattering at blue wave¬ 
lengths. ISouthworth et al.l (12015h found a strong increase in 
absorption toward the blue for the hot Jupiter WASP-103b, 
which is currently unexplained. The ac hieved precision o f the 
broad-band sp e ctropho tometry papers of iBento et alJ (12014l) and 
iMancini et al.l (l2013al) did not yield significant results for the 
hot Jupiters WASP-17b and WASP-19b, but tentatively it was 
possible to favor a certain atmospheric model before the oth¬ 
ers. In contrast, there is a larger number of broad-band spec¬ 
trophotometry investigations of hot Jupiters that did not reach 
the necessary precision to di s tinguish between atmospheric mod- 


els (ISou 

thworth et al. 20121; Nikolov et al.l 120131; Manci 

ni et alJ 

201 3bl 

2014bllal; 

Connerwheat et al. 20131; IChen et al. 

2oi4 

Eukui et al.ll20l4 

Zhou et al.||2014). 


In this work, we inve stigate the hot Jupiter HAT-P-12b. It 
was discovered in 2009 bv iHartman et al.l (l2009l) and has a mass 
of ~ 0.2 M j and a radius of about 1 Ry. The host star HAT-P-12 
is a K5 main sequence star with a mass of 0.73 Mp, a radius of 
0.7 R(7), and a magnitu de of V = 12.8 mag (iHartman et al.ll200^ 
lEhrenreich & Des^l201 Ih . The planetary mean density and the 
surface gravitational acceleration is rather low (g ~ 6 ms“^). 
This low surface gravity and the equilibrium temperature of 
~ 1000 K result in an atmospheric pressure scale height of about 
600 km. Together with the large transit depth due to the small 
stellar size, HAT-P-12b is among the known exoplanets with the 
strongest potential transmission signal. The planet orbits its host 
star ev ery ~3.2 days. HAT -P-12 was described as moderately 
active (iKnutson et aLlBOlOl) . However, the jitter of the RV mea¬ 
surem ents is low with a few meters per second (IHartman et al.l 
l2009l) . The line cores of the Ca H&K lines show emission fea¬ 
tures with log(Rj^^) =-5.1. 

Eollow - up ph otometry of HAT-P-12b was published by 
ISada et akl (12012l) . who present one J band transit light curve. 
Shortly after, t hree optical light curves were analyzed by 
iLee et all (12012l) . allowing f or a transit parameter refinement and 
a more accurate eDhemeris. lTodorov et al.l (l2013li published up¬ 
per limits on the eclipse depth and the brightness temperature at 
3.6 and 4.5 Atm. A first transmission spectrum was presented by 
iLine et all (120131) using HST in the NIR. 

The large number of transit light curves analyzed in this 
work provide the opportunity to analyze the data for transit tim¬ 
ing variation (TTV). This technique allows the detection of so 
far unseen planetary mass companions and derivation of vari¬ 
ous system parameters (lAgol et al.ll2005l) . In general, the current 
understanding of planet formation and migration predicts that 
hot Ju piters orbit their host stars alone (see, e.g., ISteffen et al.l 
12012l and references therein). Howev er, there might be e xcep- 
tions from this general rule, as when ISzabo et"^ (12013l) pre¬ 
sented several exoplanet candidates detecte d by TTVs of hot 
Jupiter s bas ed on Kepler da t a. Furt hermore. iMacieiewski et al] 
(12013l) and Ivon Essen et al.l (12013l) have announced tentative 


detections of TTV signals using ground-based data. However, 
ground-based results should be used with care as the debate 
for TTVs in the WASP-10 system illustrates (IMacieiewski et all 
1201 ll: iBarros et al.l l2013l) . There are also a handful of exam¬ 
ples of candidates or confirmed multiplanet systems with a hot 
Jupiter, where indications or evidence of anot her companion 
was ac hieved with the radial v elocity method, s ee lHartman et al.l 
(12014l) and references therein. iLee et al] (12012l) found variations 
in the transit times of HAT-P-12b that are significantly larger 
than their error values. We use our sample of transit observations 
to follow up on this potential TTV detection. 

Transit light curves always suffer to some degree from the 
correlated noise caused by t he Earth’s atmosphere or instrumen¬ 
tal effects (iPont et al.ll2006l) . One approach to mitigating this is 
to observe transits at multiple epochs and multiple telescopes 
since these light curves are not correlated to each ot her, and the 
influe nce of the individual red noise averages out (iLendl et al.l 
l2013l) . In this work we combine a large number of transits in 
their majority taken with small-to-medium class telescopes at 
different epochs to refine the ephemeris and orbital parameters 
and to search for transit timing variations and a wavelength de¬ 
pendence of the planetary radius. Sections 2 and 3 describe the 
observations and the data reduction. Section 4 presents the re¬ 
sults of the re-analysis of the literature data. The analysis and re¬ 
sults of the new transit light curves are given in Section 5, while 
Section 6 details the results of the monitoring program. Section 7 
discusses the results, followed by the conclusions in Section 8. 

2. Observations 

We acquired data of 20 primary transit events. Observations were 
done with the STELLA telescope, the Large Binocular Tele¬ 
scope (LBT), the William Herschel Telescope (WHT), the Tele- 
scopio Nazionale Galileo (TNG), the Nordic Optical Telescope 
(NOT), the Asiago 1.82m telescope, and the 0.8m telescope of 
the Montsec Astronomical Observatory (OAdM). The transits of 
the LBT and the WHT were observed in multichannel mode. The 
STELLA telescope has been used to monitor the host star HAT- 
P-12 in two colors for a length of about four months in 2012 and 
four months in 2014. A summary of the transit observations is 
given in Table [T] including the date of observation, filter, expo¬ 
sure time, time between consecutive exposures (cadence), num¬ 
ber of data points, the rms, the /J factor (see Section |4]i, and the 
airmass range during the observation. 

2.1. Transit observations 

2.1.1. STELLA 

STELLA observed 13 transits in the observing seasons 2011, 
2012, and 2014. STELLA is a robotic L2m twin telescope lo¬ 
cated on Tenerife (IStrassmeier et al.ll2004|). The instrume nt of 
use was its wide-field-imager WiESIP dWeber et al.ll2012l) with 
a field of view (EoV) of 22’x22’ on a scale of 0.32"/pixel. The 
detector is a single 4096x4096 back-illuminated thinned CCD 
with 15yL(m pixels. The telescope was slightly defocused during 
the observations (FWHM ty pically 2-3 arcseconds) to minimize 
flat-fielding errors (see, e.g.. lSouthworth et al.ll2009l) . The obser¬ 
vations have been carried out using the Sloan r’ filter in the 2011 
season and quasi-simultaneously with alternating filters Johnson 
B and V in 2012. The B band light curves are of low quality be¬ 
cause of the intrinsic faintness of HAT-P-12 at blue wavelengths 
and weather problems, therefore they are not analyzed here. In 
the year 2014, two transits were observed in Sloan g’ and one in 
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Sloan i’. The r’ observations from 2011 were performed reading 
out the entire detector area. In the season 2012, we optimized 
the setup with the read-out of a central window of Ik by Ik to 
save read-out time. Typically, the object centroid drifted by a few 
pixels per hour over the detector. 

2.1.2. LBT 

One tr ansit was observed w ith the Large Binocular Camera 
(LBC, iGiallongo et Zll2008h on April 4, 2013 (PI: I. Pagano). 
The LBC consists of two prime-focus, wide-held imagers 
mounted on the left and right arms of the LBT, and is optimized 
for blue and red optical wavelengths, respectively. We applied 
the Vspec filter on the blue arm, a filter with Sloan u’ response 
but having an increased efficiency, centered at Ac = 357.5 nm, 
and the F972N20 filter on the red arm, which is an intermediate- 
band filter centered at Ac - 963.5 nm. The two filters have been 
chosen to cover a very wide wavelength range; furthermore, the 
F972N20 filter avoids most of the telluric lines in the red region. 
On both sides, just the central chip of four chips has been read 
out to save hard disk space and a region of interest (window) has 
been defined to save read-out time. The final FoV encompassed 
~ 7.5x7.5 arcminutes^. The blue arm was defocused to an arti¬ 
ficial FWHM of about four arcseconds, the PSF of the red side 
was donut-shaped with a width of about eight arcseconds at half 
the peak count rate. The exposure time was set to 90 seconds 
on both sides with a cadence of ~ 125 seconds on the blue and 
~ 117 seconds on the red side. Unfortunately, the observation 
was affected by clouds, especially in the second half of the tran¬ 
sit, almost no flux of the target was received. However, the obser¬ 
vation was not interrupted, and a number of useful frames were 
again recorded after egress. The telescope was passively track¬ 
ing during the time series, resulting in a centroid drift of about 
20 pixels in both directions on the blue side and 80/40 pixels in 
x/y on the red side. 

2.1.3. WHT 

We observed one transit on March 15, 2014 (PI: V. Nascim- 
beni) with the triple beam, fra me-transfer CCD camera UL¬ 
TRACAM (iDhillon et alJ 1200'^ . The instrument splits the in¬ 
coming light into three channels with individual CCDs, allow¬ 
ing simultaneous observations with three filters. We chose the 
Sloan filters u’, g’, and r’. The exposure time for all three chan¬ 
nels was 5.7 seconds with almost no overheads (25 ms) owing 
to the frame-transfer technique. The telescope has been defo¬ 
cused to a FWHM of approximately 5/4/3.5 arcseconds in the 
u’/g’/r’ band. The observing conditions were photometric. The 
object centroids moved on the CCD over the time series by about 
five pixels in u’ and less than two pixels in g’ and r’. ULTRA¬ 
CAM shows no measurable flexure, therefore we attribute the 
centroid drift to atmospheric refraction. No windowing was ap¬ 
plied to the FoV of 5x5 arcminutes^. For the u’ channel, four 
exposures were co-added on chip before read-out to reduce the 
influence of the read-out noise. 

2.1.4. TNG 

One transit was observed with the TNG on May 28, 2014 using 
the focal reducer DOLORES in imaging mode (DDT program, 
PI: M. Mallonn). DOLORES is equipped with a 2048x2048 
E2V 4240 CCD camera optimized for blue wavelengths and 
a EoV of 8.6x8.6 arcminutes^. We employed the permanently 


mounted Sloan u’ filter. A binning of 2x2 pixels was used to 
reduce the read-out time. An exposure time of 40 seconds re¬ 
sulted in a cadence of 49 seconds. A mild defocus was applied 
throughout the time series to reach an artificial FWHM of about 
two arcseconds. The observing conditions were non-photometric 
with variable seeing and the target flux varied by up to 30% on 
a short time scale. Twice during the time series, the position of 
the target jumped abruptly by about ten pixels. In between these 
jumps the object drifted by less than three pixels/hour. 

2.1.5. NOT 

One transit was observed with the NOT on March 15, 2014 as 
a Fast-Track program (PI: M. Mallonn). The imaging time se¬ 
ries was taken with the stand-by CCD imager StanCam, which 
hosts a Ikxlk TK1024A CCD detector. The FoV amounted to 
3x3 arcminutes^, sufficient to image the target and one nearby 
reference star of similar brightness. We used a Bessel B filter 
with an exposure time of 60 seconds, resulting in a cadence 
of 103 seconds. The night was photometric with a seeing be¬ 
tween 1.0 and 1.3 arcseconds. We noticed an object drift of about 
six pixels in both directions over the length of the time series. 

2.1.6. Asiago 

Two transits were observed with the Asiago 1.82m telescope (PI: 
G. Piotto) and its Asiago Faint Object Spectrograph and Camera 
(AFOSC). The camera hosts a 2kx2k E2V 42-20 CCD, which 
images the 9x9 arcminutes^ FoV with a 0.26”/pixel plate scale. 
We applied a 4x4 pixel binning and read out only a region of 
interest to reduce the overheads to about two seconds. All im¬ 
ages were exposed for ten seconds. A defocus was used, which 
enlarged the PSF to a three to five arcsecond FWHM. 

2.1.7. OAdM 

One partial transit was observed on May 15, 2014 (PI: E. Her- 
rero) with the fully robotic 80-cm Ritche y-Chretien telescop e 
of the Montsec Astronomical Observatory dColome et alj|200^ . 
The telescope hosts a FLI PL4240 2kx2k camera with a plate 
scale of 0.36 arcseconds per pixel. Because of the intrinsic faint¬ 
ness of the K dwarf HAT-P-12 at blue wavelengths, we chose 
the reddest filter available here, the Cousin I filter. The exposure 
time amounted to 120 seconds with a cadence of 136 seconds. 
The beginning of the astronomical twilight prevented us from 
observing a full transit. The telescope was mildly defocused to 
an FWHM of ~ 2.5 arcseconds, the centroids drifted over the 
time series by 60 and 30 pixels in the x and y directions. 

2.2. Observations of the monitoring program 

We monitored HAT-P-12 with STELLA/WiFSIP in 2012 and in 
2014. In 2012, the observations spanned 112 from April to Au¬ 
gust. If weather, target visibility, and scheduling constraints al¬ 
lowed for it, a sequence of seven exposures in Johnson V and 
seven exposures in Johnson I was obtained every night. On a 
few occasions, the robotic telescope system observed multiple 
1+1 blocks per night separated by a few hours. In total, we gath¬ 
ered 636 images in V and 631 images in I. The exposure time 
was 30 and 60 seconds, respectively. 

In 2014, the monitoring program was continued from Febru¬ 
ary until July. At that point, 1081 images were taken in V, and 
1077 images in I. The exposure time was 60 seconds in V and 
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Fig. 1. Light curves of Table[T]after detrending with a second-order polynomial over time, presented in a chronological order according to their 
observing date. The residuals are shown in the right panel. 
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Table 1. Overview about new transit observations of HAT-P-12b. The columns give the observing date, the used telescope, the chosen filter, the 
exposure time, the observing cadence, the number of observed individual data points, the dispersion of the data points as root-mean-square (rms) of 
the observations after subtracting a transit model and a detrending function, the p factor (see Section|4]l, and the airmass range of the observation. 


Date 

Telescope 

Lilter 

^exp (s) 

Cadence 

(s) 

-^data 

rms 

(mmag) 

P 

Airmass 

2011 Mar 9 

Asiago 

R 

10 

12.1 

1117 

2.00 

1.12 

1.00-1.20 

2011 Apr? 

STELLA 

r’ 

40 

78 

205 

1.92 

1.64 

1.03-1.91 

2011 Apr 23 

STELLA 

r’ 

40 

78 

85 

1.45 

1.00 

1.07-1.28 

2011 May 9 

STELLA 

r’ 

40 

78 

204 

1.26 

1.02 

1.03-1.59 

2011 May 25 

STELLA 

r’ 

40 

78 

204 

1.83 

1.00 

1.08-2.77 

2011Jun7 

STELLA 

r’ 

40 

78 

205 

1.14 

1.00 

1.03-1.61 

2011Jun 23 

STELLA 

r’ 

40 

78 

203 

2.03 

1.06 

1.07-2.76 

2012 Eeb 6 

STELLA 

V 

30 

145 

133 

4.10 

1.08 

1.03-1.90 

2012 Mar 19 

Asiago 

R 

10 

12.1 

930 

2.11 

1.11 

1.03-1.42 

2012 May 3 

STELLA 

V 

40 

186 

92 

2.04 

1.03 

1.03-1.36 

2012 May 6 

STELLA 

V 

40 

186 

107 

2.77 

1.05 

1.03-2.02 

2012 June 4 

STELLA 

V 

40 

128 

129 

3.88 

1.00 

1.03-2.12 

2013 Apr 5 

LBT 

^ spec 

90 

125 

66 

1.57 

1.18 

1.02-1.28 

2013 Apr 5 

LBT 

L972N20 90 

117 

63 

0.98 

1.00 

1.02-1.28 

2014 Eeb 27 

STELLA 

g’ 

80 

104 

140 

1.69 

1.00 

1.07-2.17 

2014 Mar 15 

STELLA 

g’ 

80 

104 

136 

1.34 

1.45 

1.03-1.28 

2014 Mar 15 

NOT 

B 

60 

103 

139 

1.52 

1.00 

1.04-1.34 

2014 Mar 15 

WHT 

u’ 

23 

23 

676 

7.16 

1.36 

1.03-1.38 

2014 Mar 15 

WHT 

g’ 

5.7 

5.7 

2717 

1.69 

1.96 

1.03-1.38 

2014 Mar 15 

WHT 

r’ 

5.7 

5.7 

2717 

1.27 

2.00 

1.03-1.38 

2014 May 15 

OAdM 

I 

120 

136 

40 

1.44 

1.76 

1.20-1.51 

2014 May 28 

STELLA 

i’ 

45 

72 

240 

1.55 

1.00 

1.04-1.56 

2014 May 28 

TNG 

u’ 

40 

49 

308 

5.44 

1.22 

1.03-1.32 


50 seconds in I. The images were observed in blocks of five ex¬ 
posures per filter. In the year 2013, the CCD detector of WiF- 
SIP was replaced, therefore the time series in 2012 and 2014 
were not obtained with identical instruments. Additional tele¬ 
scope manufacturing work was done in May 2014, followed by 
the development of a new pointing model, which manifests in 
the time series as an abrupt change in detector position on June 

2. Unfortunately, this change was accompanied by a jump in the 
differential light curve of about one percent amplitude. 

3. Data reduction 

3.1. Transit data 

The reductions of the transit data obtained on STELLA, LBT, 
WHT, TNG, NOT, and OAdM were done homogeneously with 
a customized ESO-MID AS pipeline, which call s the photome¬ 
try software SExtractor (iBertin & Arnoutslll996l) . Bias and flat- 
field correction was done in the usual way, with the bias value 
extracted from the overscan regions. Cosmic ray correction was 
carried out within SExtractor. We performed aperture photome¬ 
try with circular apertures (MAG_APER in SExtractor) and au¬ 
tomatically adjusted elliptical apertures (MAG_AUTO in SEx¬ 
tractor). The elliptical aperture offers the advantage of account¬ 
ing for variations in the PSE over time, but with the disadvantage 
of different apertures for target and comparison stars. 

Our algorithm automatically searches for the combination of 
comparison stars, which minimizes the rms of the target light 
curve. The comparison stars are weighted according to their indi- 
vidual light curve quality in terms of point-to-point scatter (rms) 
(iBroeg et al.l2005h . The rms of the target light curve is computed 


not only on the out-of-transit part, but also on the residuals of 
the entire time series after subtracting a transit fit including a 
second-order polynomial in time for detrending (see SectionUJi. 
As expected, the algorithm tends to choose only a low number of 
nearby stars in non-photometric conditions because atmospheric 
extinction differs for stars spatially separated on the sky. In con¬ 
trast, in photometric conditions, the number of chosen compari¬ 
son stars is higher (typically 5 to 8), including stars with larger 
spatial distance. 

In a last step, the algorithm repeats the light curve extrac¬ 
tion many times with different aperture sizes to search for the 
aperture that minimizes the scatter in the differential target light 
curve. In the majority of our observations, the circular aperture 
with fixed radius (MAG_APER) yielded the more stable pho¬ 
tometry than the flexible elliptical aperture (MAG_AUTO). 

The LBT transit observation suffered from highly variable 
sky transparency due to clouds moving through the EoV. We 
restrict the time series to data points that reached at least half 
the maximum flux value. We carried out a fringe correction for 
the LBC E972 light curve that resulted in a marginal depression 
of the point-to-point scatter of the final transit light curve from 
1.01 mmag to 0.98 mmag. The values of the transit parameters 
remained unchanged. 

The Asiago data were redu ced with the software pa ckage 
STARSKY, described in detail in iNascimbeni et al.l(l2013a|l . Eig- 
ure[T] shows the differential light curves of all transits after nor¬ 
malization and detrending with a second-order polynomial over 
time. 
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3.2. Monitoring data 

For the monitoring data, the bias and flatfield correction was 
done with the STELLA data reduction pipeline (for details 
see Granzer et al. 2015, in preparation). The following reduc¬ 
tion steps were done with routines written in ESO-MIDAS. We 
conducted aperture photometry with SExtractor applying the 
MAG_AUTO option to automatically adjust elliptic apertures. 
This method provides the flexibility to account for the varying 
observing conditions over the ten months of observing time. The 
I band data were fringe-corrected. We verified that the three used 
comparison stars GSC2.2-N130301273, GSC2.2-N130301284, 
and GSC2.2-N130301294 were photometrically stable to our 
achieved precision. 

We discarded data points of very high airmass, low flux 
caused by poor weather, extremely high sky background (twi¬ 
light), extreme values of FWHM indicating problems with the 
telescope focus, and extreme PSE elongation that indicates vi¬ 
brations of the telescope caused by wind. Eurthermore, we ex¬ 
cluded the data taken after June 2, 2014 with the new pointing 
model, as the new detector position caused a jump in the differ¬ 
ential light curve of about 1.0% in V and 0.5% in I. Additionally, 
we applied a 5 cr clipping to remove persistent outliers. The final 
light curve in V contains 486 of the 636 observed data points in 
2012 and 873 of the 1081 data points in 2014. The I band light 
curve contains 402 out of 631 data points for 2012 and 709 out 
of 1077 data points for 2014. 

Einally, the V band light curve of 2012 has a length of 
109 days with a scatter of 5.1 mmag (rms), while the data points 
of 2014 cover 112 days and show a deviation of 3.2 mmag. The 
deviation is lower simply because the exposure time was doubled 
in 2014. The I band light curve from 2012 scatters by 3.0 mmag, 
which is comparable to the scatter of 2.6 mmag of the I band 
data from 2014. 


4. Transit modeling and results 

4.1. Giobai transit parameters 


Throughout this work, we model the transit light curves with 
the p ublicly available software JKTEBOFQ (ISouthworth et al.l 
I 2 OO 4 I 1 . It allows for a simultaneous fit for the transit model and 
detrending. Most ground-based differential light curves show a 
long-term trend potentially caused by atmospheric effects due to 
a color difference between the target and the combined reference 
stars. Throughout this HAT-P-12b analysis, we use a second- 
order polynomial over time to detrend the light curves. The tran¬ 
sit fit parameters consist of the sum of the fractional planetary 
and stellar radius, r* -H and their ratio k - rpjr./^, the or¬ 
bital inclination i, the transit midtime Tq, and the host-star limb- 
darkening coefficients (LDC) u and v of the quadratic limb dark¬ 
ening law. The index refers to the host star and “p” refers 
to the planet. The dimensionless fractional radius is the absolute 
radius in units of the orbital semi-major axis a, r* = R*/a, and 


= R 


nia. The plane tary eccentricity is fixed to zero following 


iHartman et ah (l2009h a nd the orbital period Porb to 3.21306 days 
according to Lee et al.l (l2012h . Additional fit parameters are the 
three coefficients co,i ,2 of the parabola over time. 

Throughout this work, we model the stellar limb darkening 
with a quadratic limb darkening law. Initial values for the LDC 
u and V are found by inter polating tabulated ATLAS values from 
IClaret & Bloemeril (1201 111 for the stellar parameters of HAT-P- 
12: T = 4650 K, log g = 4.61, [Pe/H]= -0.29, and v„„ 2 . = 


* http://www.astro.keele.ac.uk/jkt/codes/jktebop.html 


0.85 km s ' (IHartman et al.l20()^ . In the transit modeling, we fit 
for the linear coefficient u but fix v to its theoretical value. This 
ensures enough flexibility that the transit fits are not significantly 
biased by systematic errors of the LDC calculation. On the other 
hand, it avoids parameter degeneracies that occur when fitting 
for both u and v. We account for errors in the calculated v by a 
perturbation of this value by +0.1 during the error estimation of 
the transit parameters. 

Eor the error estimation, we run a Monte Carlo simula¬ 
tion (“task 8” in JKTEBOP) o f 5000 steps and the residual- 
perm utation algorithm (“task 9”. lJenkins et al.l2002l : ISouthworthl 
l2008h and choose the higher value. The latter method takes the 
presence of correlated noise into account. 


4.1.1. Fit of indivicdual transits 

In a first step we apply a 4 cr rejection to the light curve resid¬ 
uals after subtracting an initial model, with cr as the standard 
deviation. A subsequent fit is used to rescale the individual pho¬ 
tometric error bars: The error bars supplied by SExtractor yield 
a reduced that is always slightly greater than 1, which indi¬ 
cates that the photometric error bars are underestimated. There¬ 
fore, we enlarge the error bars by a factor individual to each light 
curve to give a of unity. Eurtherm ore, we calcu l ate th e so- 
called 13 factor, a c oncept introduced bv iGillon et al.l (l2006h and 
IWinn et ^ (l2008l) to include the contribution of correlated noise 
in the light curve analysis. It describes the evolution of the stan¬ 
dard deviation cr of the light curve residuals when they become 
binned in comparison to Poisson noise. In the presence of cor¬ 
related noise, ctn of the binned residuals is larger by the factor 
ff than with pure uncorrelated (white) noise. There is a depen¬ 
dence of this value on the bin size, and generally the duration 
of ingress/egr ess is used as a re levant time scale for transit pho¬ 
tometry (e.g.. lWinn et al.ll2008h . We estimate the p factor given 
in Tables [1] and 13 as the average of the values from 10- to 25- 
minute bins, since the ingress of HAT-P-12b lasts about 19 min¬ 
utes. It may happen that the derived p factor is slightly smaller 
than unity because of statistical fluctuations. Such cases were 
manually set to p - \ .0. The derived p factor for each transit 
light curve is used to further enlarge the individual photometric 
error bars, and the final transit fit was performed including error 
estimation. 

We can compare the scatter of the transit parameters de¬ 
rived per light curve with their one-sigma uncertainties to draw 
a rough conclusion about the reliability of the uncertainties: The 
25 set of parameters from the complete transits analyzed in this 
work give a reduced x^ for the inclination i of 2.2, which can 
be interpreted as slightly underestimated error bars. Eor the frac¬ 
tional stellar radius r*, the reduced;^^^ is 1.0 and 0.9 for the linear 
LDC u (derived values of u were compared with their theoreti¬ 
cal predictions), an indication for error bars that match the scat¬ 
ter within these 25 measurements very well. The reduced for 
the radius ratio k is rather high with a value of 3.5, so that ei¬ 
ther the errors are underestimated, or else we see intrinsic effects 
such as a potential wavelength dependence (see Section 1431) or 
light curve modifications caused by starspots (see Section|5]l, for 
which k is the most sensitive transit parameter. 


4.1.2. Re-analysis of the literature data 

We analyze the published transit light curves described in Table 
13 again. Before we include the literature light curves in our data 
sample, we test whether our analysis routines are able to repro- 
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Table 2. Overview of published transit observations of HAT-P-12b, 
which are re-analyzed in this work. 


Date 

Filter 

^data 

rms 

(mmag) 

yS 

Reference 

2007 Mar 27 

i’ 

151 

1.78 

1.00 

Hartman et al. t2009) 

2007 Apr 25 

z’ 

375 

2.62 

1.02 

Flartman et al. (2009) 

2009 Feb 5 

z’ 

219 

2.66 

1.04 

Hartman et al. (2009) 

2009 Mar 6 

g’ 

214 

5.14 

1.31 

Hartman et al. t2009) 

2010 May 31 

j 

383 

2.52 

1.33 

Sadaetal. (2012) 

2011 Mar 29 

R 

224 

1.75 

2.24 

Lee et al. (2012) 

2011 Apr 14 

R 

273 

2.64 

1.68 

T.ee et al. [2(|)1 2') 

2011 May 13 

R 

314 

1.80 

1.57 

Lee et al. (2012) 


duce the published results. Whenever possible we treat the pub¬ 
lished data identically to ou r observations (see Section 14.1.1b . 
For the iHartman et al.l (l2009h light curves, we achieve excellent 
agreement in the paramete rs, but find error bars that are higher 
by up to a factor of 2. The lLee et al.l (1201 2l) results, which were 
also derived with JKTEBOP, are reproduced as well, however 
our error estimation yields values higher by a factor of 2 for the 
planet-to-star radius ratio and 10 for the orbital inclination. The 
source of th is discrepancy rem ains unexplained. The J-band light 
curve from ISada et ahl (1201 2h shows rather complex deviations 
from the symmetric transit shape and was originally detrended 
by a fifth-order polynomial over time fitted to the out-of-transit 
data. The uncertainty of this polynomial was not included in the 
transit parameter error estimation. Following this approach, we 
can reproduce the transit parameters to within 1 cr with about the 
same error bars. However, a simultaneous fit of detrending poly¬ 
nomial and transit model changes the parameters by up to 1.5 cr 
and increases the uncertainty by up to a factor of 2. 


4.1.3. Fit of multiple transits 

We next fit all light curves of one filter band simultaneously, 
again using the inflated photometric error bars of Section |4. 1.1 1 
JKTEBOP does allow the simultaneous fit of multiple transit ob¬ 
servations if one joint limb-darkening function can be applied. 
The software cannot attribute different LDC to subsamples of 
the input data, which prohibits the simultaneous fit of multicolor 
data. However, the current version 34 of JKTEBOP is able to as¬ 
sign independent detrending polynomials to the individual light 
curves, so the detrending is included in the transit fit. Eree pa¬ 
rameters of the fits are rp + r*, k, i, Tq, the period Porb, the linear 
coefficient u of the quadratic law, and a set of coefficients co,i ,2 
of the detrending polynomial for each individual transit obser¬ 
vation. The quadratic LDC v is fixed to its theoretical value, but 
perturbed during the error estimation. Eor the B band transit, the 
only filter with a single transit measurement in Table [3] we also 
keep u to its theoretical value to ensure that the fit converges 
at reasonable parameter values. Table [3 lists the transit param¬ 
eters derived from all complete transits of each band fitted si¬ 
multaneously. Eor a simultaneous fit, it is necessary to convert 
all barycentric Julian Date (BJD) timings of the observations 
to a common time standard which is in this work the barycen- 
tric dynamica l time (TDB) following the recommendation of 
[Eastman et alJ (l2010h . The u’ band transits are not listed because 
their rather low accuracy would not improve the final results. The 
final transit parameter values are formed as the weighted average 
of the seven filters. 


Eor all our filter bands with multiple transits, the derived un¬ 
certainties of the simultaneous fit are larger and more conser¬ 
vative than the uncertainties derived by a weighted average of 
the individually fitted transits (cr^ = l/S"^jCr“^). A plausible 
explanation is that the error estimation becomes more sensitive 
to correlated noise when analyzing multiple transits simultane¬ 
ously. This is indicated by the fact that the residual-permutation 
method gives larger uncertainties than the Monte Carlo simu¬ 
lation. However, the uncertainties in the seven bands seem too 
conservative in part; the final values for i and r* yield a of 
0.2 and 0.1, respectively. Eor u, we derive s^Xied 
the x^ed ^ decreases from 3.5 for the 25 individual complete 
transits to 1.3 for the seven per-band values. 

Our final transit parameters are generally in good agreement 
with previous publications on HAT-P-12b. Most values agree 
within the 1 cr levels, and the strongest deviatio n of about 2cr 
is foun d for k between our value and t he one from Hartman et al.l 
(I 2 OO 9 I 1 . The rather high value for i of iLee et"^ (1201 2h matches 
the value of this work to 1 cr if we consider the uncertainty 
of +1.0 deg, which we derive from a re-analysis of their light 
curves. 


4.2. Ephemeris and 0-C diagram 

To derive the ephemeris and investigate potential transit timing 
variations, we fit all individual transits, complete or partial, with 
a transit model fixing + r*, k, and i to the final parameter 
values of Table [3] Eree parameters are the transit midtime To, 
the linear LDC u and the three coefficients of the detrending 
polynomial co.i.i- The uncertainties are estimated using a Monte 
Carlo simulation with 500 0 steps. Although the light curves from 
IHartman et al.l (l2009l) and iLee et al.1 (120121) used slightly differ¬ 
ent detrending functions, we approximate their error contribu¬ 
tion by including the same three free parameters co,i ,2 as for the 
newly observed light curves of this work. The fit and MC error 
estimation of the Sada J band light curve additionally contain 
the free parameters C3, C4, and C5 to model the trend with a fifth- 
order polynomial over time. The quadratic LDC v is fixed to its 
theoretical value. We do not perturb v in the MC simulation since 
we find that To is insensitive to the treatment of limb darkening. 

Table |4] summarizes all individual transit timings with their 
uncertainties and includes the individual deviation O - C of the 
observation from the calculated ephemeris. All times are given 
in BJD(TDB). We use a linear least-squares fit to the transit 
times of the 26 complete transits to newly determine a linear 
ephemeris: 


r, = BJD(TDB)2456032.151332(39) + 3.21305756(20) A, 

( 1 ) 

where Tc is the predicted central time of a transit, N the cycle 
number with respect to the reference midtime, and the num¬ 
bers in brackets give the uncertainties of the last two digits es¬ 
timated with a Monte Carlo simulation of 5000 steps. The ref¬ 
erence midtime is chosen to minimize the covariance between 
reference midtime and period. The estimated orbital period dif¬ 
fers b y about 1 cr from the original value given in the discovery 
paper IHartman et ^ (l2009l) . but is an order of magnitude more 
accurate. Howe ver, our p e riod i s 0.15 seconds shorter than the 
period value of iLee et"^ (1201 2l) . which is a deviation of about 
5 cr. We note that some of the very first transits of the data sam¬ 
ple, the three complete Hartman transits, happ en systematically 
earlier by about 65 seconds in the calculation of iLee et^(l2012h 
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Table 3. Transit parameters of HAT-P-12b derived in this work and comparison to previously published values. The final values were formed as 
the weighted average of the seven filters. 


Filter 

r* = R*/a 

k = rp/r* 

i(deg) 

linear LDC u 

quad. LDC v 

B 

0.08526 + 0.00382 

0.13916 + 0.00311 

89.00+ 1.04 

0.9150" 

-0.0592 

g’ 

0.08584 + 0.00251 

0.13847 + 0.00215 

88.60 + 0.48 

0.8080 + 0.0371 

-0.0122 

V 

0.08510 + 0.00602 

0.13679 + 0.00439 

89.41 + 1.24 

0.7458 + 0.0972 

0.0961 

r’ 

0.08456 + 0.00145 

0.13611 + 0.00150 

89.06 + 0.78 

0.5829 + 0.0299 

0.1665 

R 

0.08565 + 0.00155 

0.13826 + 0.00167 

89.07 + 0.66 

0.5790 + 0.0377 

0.1823 

i’ 

0.08551 +0.00282 

0.13515 + 0.00230 

89.26+ 1.01 

0.5081 +0.0603 

0.2198 

z’ 

0.08561 + 0.00290 

0.14082 + 0.00204 

88.92 + 0.88 

0.3079 + 0.0590 

0.2426 

Final 

0.08532 + 0.00085 

0.13779 + 0.00079 

88.98 + 0.29 



Hartman et al. 120091 

0.0850 + 0.0013 

0.1406 + 0.0013 

89.0 + 0.4 



Lee et al. 120121 

0.0852 + 0.0012 

0.1370 + 0.0019 

89.9 + 0.1 



Sada et al. 120121 

0.0891 + 0.0045 

0.1404 + 0.0026 

88.5 + 1.0 



Line et al. 120131 

0.0862 + 0.0029 

0.137 + 0.0011 

88.7 + 0.6 




Notes. Both limb darkening coefficients are fixed to its theoretical values for filter B, see text for details. 


than in our calculation. Since the data of iHartman et al.l (l2009h 
are given in the time standard UTC, this deviation might be ex¬ 
plained by a faulty UTC-TDB correction. 

Simultaneously observed transits offer the possibility of test¬ 
ing the reliability of the derived transit timing uncertainties, al¬ 
though one needs to be careful to draw definite conclusions from 
such low number statistics. We observed the transit of March 15, 
2014 with three different telescopes, the William Herschel Tele¬ 
scope observed with ULTRACAM in three different bands. We 
therefore have five transit light curves of the very same event. 
The reduced value of these five transit midtimes is 2.32. We 
ran a simulation of 10 000 datasets of five points each, added 
Gaussian noise, and find that such a higher is reached 

in only 10% of all cases. Therefore, it appears likely that the 
uncertainties are slightly underestimated. However, with a hy¬ 
pothetical increase of only 10% on the error bars, thex^ed 
would already lie within 1 cr of the mean of the simulated x^^^ 
distribution. 

Figure |2] shows the deviations of the individual transit times 
from the linear ephemeris of Equation [1] All data points give 
X^^^ ~ 4.8, and an exclusion of the partial transit lowers this 
value only marginally to;^^^^J ~ 4.5. The most deviating point is 
seven sigma away from the predicted ephemeris and corresponds 
to the March 19, 2012 transit from the Asiago telescope. If this 
point is excluded, the^^^^ lowers to 2.4. 

We calculate a Lomb-Scargle periodogram of all transits 
times excluding the partial transits shown in Figure[3] For a test, 
we again exclude the Asiago transit from March 19, 2012 and 
find a very different periodogram without any significant peak. 
The limit for the highest frequency is defined by the average 
Nyquist frequency. We conclude that we find no evidence of 
any transit timing variation and that the rather high value of 
X^^^ ~ 2.4 is caused by an underestimation of the error bars 
due to systematics. One potential source of syste matics in the 
light curves that influences the timing are starspots (iBarros et al.l 
1201 3h . a plausible scenario for this slightly variable host star (see 
Section]?]). 


Table 4. Observed transit times of HAT-P-12b. 


BJD(TDB) 

(2,450,000+) 

Uncertainty 

Cycle 

number 

0-C 

Reference 

4187.85647“ 

0.00061 

-574 

0.00017 

Hallman et al. 120091 

4216.77343“ 

0.00023 

-565 

-0.00037 

Hartman et al. ('2009) 

4869.02462“ 

0.00033 

-362 

0.00013 

Hartman et al. ('2009) 

4897.94297“ 

0.00057 

-353 

0.00095 

Hallman et al. 120091 

5347.76952“ 

0.00035 

-213 

-0.00054 

Sada et al. 12012) 

5630.51929 

0.00011 

-125 

0.00015 

this work (Asiago) 

5649.79770“ 

0.00035 

-119 

0.00022 

Lee et al. (2012) 

5659.43649 

0.00044 

-116 

-0.00016 

this work (STELLA) 

5665.86227“ 

0.00037 

-114 

-0.00049 

Lee et al. (2012) 

5675.50436 

0.00047 

-111 

0.00242 

this work (STELLA) 

5691.56661 

0.00022 

-106 

-0.00061 

this work (STELLA) 

5694.78089“ 

0.00029 

-105 

0.00060 

Lee et al. (2012) 

5707.63244 

0.00027 

-101 

-0.00006 

this work (STELLA) 

5720.48443 

0.00015 

-97 

-0.00031 

this work (STELLA) 

5736.54999 

0.00030 

-92 

-0.00004 

this work (STELLA) 

5964.67787 

0.00059 

-21 

0.00075 

this work (STELLA) 

6006.44779 

0.00012 

-8 

0.00091 

this work (Asiago) 

6051.42937 

0.00048 

6 

-0.00030 

this work (STELLA) 

6054.64202 

0.00060 

7 

-0.00070 

this work (STELLA) 

6083.55936 

0.00065 

16 

-0.00088 

this work (STELLA) 

6388.80101 

0.00081 

111 

0.00029 

this work (LBT) 

6388.80064 

0.00048 

111 

-0.00007 

this work (LBT) 

6716.53165 

0.00033 

213 

-0.00093 

this work (STELLA) 

6732.59827 

0.00063 

218 

0.00039 

this work (STELLA) 

6732.59765 

0.00011 

218 

-0.00022 

this work (NOT) 

6732.59811 

0.00028 

218 

0.00023 

this work (WHT) 

6732.59741 

0.00023 

218 

-0.00046 

this work (WHT) 

6732.59796 

0.00009 

218 

0.00008 

this work (WHT) 

6793.64833 

0.00092 

237 

0.00236 

this work (OAdM) 

6806.49642 

0.00062 

241 

-0.00178 

this work (STELLA) 

6806.49815 

0.00019 

241 

-0.00004 

this work (TNG) 


Notes. Transit midtime newly estimated in this work. 


4.3. Planetary radius over wavelength 

One major goal of this work is to search for a wavelength depen¬ 
dence for the effective planetary radius. For this purpose, we fit 
all light curves again individually, this time by fixing the incli¬ 
nation i and the fractional stellar radius r* to the values derived 
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Fig. 2. Transit timing residuals versus the linear ephemeris of Equation 
[T] The filled symbols represent complete transits, and the open squares 
indicate transits which were just partially observed. 



frequency [l/day] 


Fig. 4. Planet-to-star radius ratio k for the individual transit light curves. 
Black crosses show k when the linear LDC n is a free-fit parameter, and 
red points show the same value when u is fixed to its theoretical value 
given in Table [5] The error bars of the latter are on average 20 percent 
smaller than the uncertainties of the former so are not shown here. Data 
points in green show the value of k for the simultaneous fit of multiple 
transits per filter. The horizontal line shows the globally derived value 
of k as given in Tablej^and its uncertainty in dashed lines. 


We combine the data sets per filter band again in a simul¬ 
taneous fit. Free parameters are m, and a set of coefficients 
co,i ,2 for every individual transit observation. The two partial 
STELLA r’ transits and the partial OAdM transit are not ana¬ 
lyzed here. In general, we note very good agreement of the ob- 
servationally derived linear LDC per filter and their theoretical 
calculations. Only in the two filter sets mostly affected by red 
noise, R and J, do we see deviations of 2.1 and 2.3 cr. 

The results for individual transits and per-filter fits are given 
together in Tableland Ligure|4] 


Fig. 3. Upper panel: Lomb-Scargle periodogram of the timing devia¬ 
tion from the linear ephemeris (O-C) when excluding the partial transits 
(black). A marginal signal can be seen at a frequency of ~0.0019 day“*. 
In comparison, the periodogram without the transit of March 19, 2012 
(epoch=-8) is also shown (red). Lower panel: Window function of the 
data set. 

in Section |4] and the orbital period and transit midtime to val¬ 
ues derived from the ephemeris (Eq.[T]i in Section IA21 All these 
parameters are assumed to be wavelength independent, but their 
values have uncertainties. Here, we neglect Ar*, A/, and AT^ un¬ 
der the assumption that they are a source of un certainty com¬ 
mon to all band passes following, for example, iMancini et alJ 
(l2014ah . The derived error bars /Sk are then relative errors. Lree 
fit parameters of the model are the fractional planetary radius 
and the three coefficients of the detrending parabola co,i, 2 - The 
J-band light curve needs the additional fit parameters C 3 _ 4_5 be¬ 
cause of its more complex trend. Lor comparison we fit each 
transit twice, first with the linear stellar LDC u as fit parameter 
and second with u fixed to its theoretical value. The quadratic 
stellar LDC v is fixed to its theoretical value in both cases but 
perturbed in the error estimations. The partial transit of the LBC 
covers too little orbital phase to reasonably constrain detrending 
and limb darkening in the same fit. Therefore, for the u’ band 
light curve, u is fixed to the value found in a simultaneous TNG 
-I- WHT u’ band fit, and for the L972N20 filter, u is fixed to the 
value derived by the z’ band light curves. 


5. Results of the monitoring program and starspot 
correction 

We monitored HAT-P-12 in 2012 and 2014 for about four 
months each in two colors, V and I, to investigate photo¬ 
metric variations potentially caused by starspots. Such bright¬ 
ness inhomogeneities on the stellar photosphere, whether oc¬ 
culted or unocculted by the transiting planet, modi fy the de¬ 
rived transit parameters, especially the tra nsit depth (ISing et al.l 
1201 Ibl^nt et al.l2013HOshagh et alj2014l) . The differential light 
curves are presented in Ligure|5] A calculated Lomb-Scargle pe¬ 
riodogram for V and I in 2012 shows no compelling periodic¬ 
ity in this dataset. The periodogram of the V band data of 2014 
shows indications of variations with a period of about 60 days. 
A fit of a sine curve with this period yields a semi-amplitude of 
1.8 + 0.4 mmag. The I-band data set of 2014 also shows mild 
indications for this period, but a period of about half this value, 
P ~ 30 days, is more significant. A sine fit yields an amplitude 
of 1.9 + 0.3 mmag for the 30-day period and 1.4 + 0.4 mmag 
for the period of 60 days. We calculated plain least-squares peri- 
odograms to estimate the amplitude of periodic sine and cosine 
functions over a wide period range and concluded that periodic 
signals of semi-amplitudes larger than about 2 mmag are ex¬ 
cluded by the V and I data of 2014. Lor the observation period of 
2012, we can exclude periodic signals of semi-amplitudes higher 
than about 1.5 mmag. 
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Table 5. Planet-to-star radius ratio k per transit observation and per filter with relative uncertainties. The parameters r^,, i, Tq, and Po* are fixed 
to their values derived in Section|4] In Column 3, k is calculated with u as free parameter and v fixed to theoretical value, but perturbed in error 
estimation. In Column 5, k is calculated with u and v fixed to theoretical values. 


Date 

Filter 

M free 

u fixed 

V 



k 

u 

k 

u 

fixed 

Individual light curves: 



2013 Mar 5 

u’ 

0.13874 ±0.00216 

1.0968“ 

0.13947 

1.0698 

-0.2218 

2014 Mar 15 

u’ 

0.13989 ±0.00282 

0.9922 ± 0.0797 

0.13793 

1.0698 

-0.2218 

2014 May 28 

u’ 

0.13627 ±0.00263 

1.1778 ±0.0741 

0.13914 

1.0698 

-0.2218 

2014 Mar 15 

B 

0.13836 ±0.00106 

0.9428 ± 0.0301 

0.13905 

0.9154 

-0.0592 

2014 Feb 27 

g’ 

0.13932 ±0.00119 

0.8046 ± 0.0346 

0.13879 

0.8492 

-0.0122 

2014 Mar 15 

g’ 

0.13328 ±0.00125 

0.9074 ± 0.0359 

0.13392 

0.8492 

-0.0122 

2014 Mar 15 

g’ 

0.13605 ±0.00049 

0.8662 ± 0.0143 

0.13620 

0.8492 

-0.0122 

2012 Feb 6 

V 

0.14103 ±0.00244 

0.7060 ± 0.0798 

0.14180 

0.6934 

0.0961 

2012 May 3 

V 

0.13615 ±0.00179 

0.6776 ± 0.0607 

0.13661 

0.6934 

0.0961 

2012 May 6 

V 

0.13821 ±0.00221 

0.8590 ± 0.0695 

0.14114 

0.6934 

0.0961 

2012Jun4 

V 

0.13783 ±0.00251 

0.6575 ± 0.0857 

0.13785 

0.6934 

0.0961 

2011 May 9 

r’ 

0.13915 ±0.00097 

0.5350 ± 0.0370 

0.13881 

0.5831 

0.1665 

2011 May 25 

r’ 

0.13951 ±0.00103 

0.6473 ± 0.0659 

0.13923 

0.5831 

0.1665 

2011Jun7 

r’ 

0.13731 ±0.00102 

0.5796 ± 0.0397 

0.13749 

0.5831 

0.1665 

2011Jun 23 

r’ 

0.14187 ±0.00122 

0.5818 ±0.0416 

0.14164 

0.5831 

0.1665 

2014 Mar 15 

r’ 

0.13377 ±0.00038 

0.6047 ± 0.0129 

0.13420 

0.5831 

0.1665 

2011 Mar 29 

R 

0.13930 ±0.00152 

0.7328 ± 0.0484 

0.14312 

0.5450 

0.1823 

2011 Apr 14 

R 

0.14176 ±0.00160 

0.6319 ±0.0579 

0.14325 

0.5450 

0.1823 

2011 May 13 

R 

0.14067 ±0.00130 

0.5322 ± 0.0767 

0.14044 

0.5450 

0.1823 

2011 Mar 9 

R 

0.13880 ±0.00064 

0.5977 ± 0.0213 

0.13984 

0.5450 

0.1823 

2012 Mar 19 

R 

0.13598 ±0.00062 

0.5806 ± 0.0161 

0.13687 

0.5450 

0.1823 

2007 Mar 27 

i’ 

0.13931 ±0.00205 

0.2295 ±0.1119 

0.14029 

0.4379 

0.2198 

2014 May 28 

i’ 

0.13597 ±0.00083 

0.4786 ±0.0310 

0.13715 

0.4379 

0.2198 

2007 Apr 25 

z’ 

0.14047 ±0.00107 

0.3261 ± 0.0381 

0.14018 

0.3518 

0.2426 

2009 Feb 5 

z’ 

0.14024 ±0.00142 

0.3633 ± 0.0532 

0.14057 

0.3518 

0.2426 

2013 Apr 5 

F972 

0.13888 ±0.00182 

0.3302* 

0.13868 

0.3518 

0.2426 

2010 May 31 

J 

0.13775 ±0.00285 

0.3660 ± 0.0637 

0.14033 

0.2200 

0.3120 

Combination of light curves per filter: 



u’ 

0.13728 ±0.00244 

1.0968 ±0.0656 


1.0698 

-0.2218 


g’ 

0.13633 ±0.00074 

0.8711 ±0.0123 


0.8492 

-0.0122 


V 

0.13857 ±0.00131 

0.7018 ± 0.0394 


0.6934 

0.0961 


r’ 

0.13741 ±0.00071 

0.5833 ±0.0130 


0.5831 

0.1665 


R 

0.13793 ±0.00098 

0.5870 ± 0.0201 


0.5450 

0.1823 


i’ 

0.13703 ±0.00103 

0.4695 ± 0.0297 


0.4379 

0.2198 


z’ 

0.14035 ±0.00087 

0.3302 ± 0.0289 


0.3518 

0.2426 


Notes. The uncertainties of k with u fixed to its theoretical value are on average 20 percent smaller than with u as free to fit the parameter, but are 
less realistic because they do not include an uncertainty of u. For the combined fit of multiple transits per filter, the fitted values of ii generally 
agree with their theoretical values. The linear limb darkening coefficient u was fixed but perturbed to the value derived from the u’ band light 
curves 2014, Mar 15 and 2014, May 28. ® The linear limb darkening coefficient u was fixed but perturbed to the value derived from the z’ band 
light curves 2007 Apr 25 and 2009 Feb 5. 


An inspection of the light curves by eye reveals a strongest 
measured variation of almost 1% in the V band data of Season 
2014 between Julian date 2456701 and 2456723. The difference 
A/ of the measured stellar flux and the flux of the unspotted 
photosphere might be larger than the measured variability, be¬ 
cause a certain level of permanently visible starspots may exist. 
We approximate this level of perman ent flux dimmin g by the 
variance of our monitoring light curve dPont et al.ll2013h and ob¬ 
tain A/ ~ 1.3%. If we assume a temperature contra st between 
photosphere and spots of 1000 K (IStrassmeieiil2009t) . the filling 
factor would amount to ~ 0.016. Correspondingly, the maximum 
correction for unocculted spots ranges from about Ak - 0.0010 


in the u’ band to about Ak - 0.0006 at J using Equation 4 in 
ISing et ail (1201 ih and the assumption of black body radiation 
for the photosphere and the spots. These correction values re¬ 
semble our 1 cr uncertainties for k in Table |5] However, this es¬ 
timation is based on the 1% measured flux variation, and none 
of our measurements took place at this moment of highest mea¬ 
sured activity. Instead, to all transits observed during the moni¬ 
toring campaigns, we could attribute filling factors of only half 
the mentioned value or less. Therefore, the contribution from an 
unocculted star spot would be significantly smaller than our un¬ 
certainties in most cases, so we do not include this correction. 
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Fig. 5. Upper panel: STELLA/WiFSIP differential V band photometry with the original data points in green and the average of the five exposures 
taken as a block given in black, with the error bars as the standard deviation divided by the square root of the number of averaged data points. 
On the left the observations from 2012 are shown, on the right the data from 2014. Middle panel: Same as in the upper panel, here for the I-band 
photometry. Lower panel: The Lomb-Scargle periodogram for the V band data (green) and the I-band data (red), on the left for 2012, on the right 
for 2014. 


Although the correction for occulted spots could in principle 
be larger, we keep the uncorrected values in our analysis for the 
following reasons: 

- At the level of precision of our light curves, we do not see 
any indication of a bump caused by a starspot crossing. 

- Transits affected mostly by occulted spots are shallower than 
the average value. The transits of 2014 March 15 show a 
shallower transit in g’ and r’. However, their simultaneously 
taken equivalents in B and u’ are not shallower than the av¬ 
erage value. Because the effect of occulted spots would be 
enhanced for these bluer wavelengths, we conclude that the 
lower values of k in g’ and r’ are caused by different kinds of 
systematics. 

- If the spots are not homogeneously distributed along the tran¬ 
sit chord, they cause deformations that is treated like red 
noise in our analysis, so that the uncertainty introduced by 
these spots is included in the uncertainty of the transit pa¬ 
rameters. In the less likely case of spots evenly distributed 
along the transit chord (resulting in a symmetric transit shape 
without bumps), the shallower transit curve would not match 
other transits taken at moments of less stellar activity and be 
considered as correlated noise in the error analysis. 

- The monitoring campaigns have shown that the activity is 
time dependent, with times of no measurable spot modula¬ 
tion. We interpret these as times with very low spot activity 
and assume that the simultaneous analysis of transits taken 
at different epochs weakens the influence of spots crossed at 
times of higher activity. 

6. Discussion 

6.1. Transmission spectrum 

We compare the measured transmission spectrum (referring to 
the common fits in each band) to three different theoretical at¬ 


mosphere models: The first model presents a cl oud-free atmo¬ 
sphere and was computed bv iFortnev et al.l (l201C)h using a planet¬ 
wide average P-T profile, solar metallicity, and the planet param¬ 
eters of HAT-P-12b. The second model simulates an atmosphere 
with a haze layer by the introduction of a Rayleigh-scattering 
slope of dRn/d\ nA - -AH similar to the observ ed spectrum 
of HD189733b (iLecavelier Des Etangs et aDl2008h . Here, H is 
the atmospheric scale height and equals ~ 600 km in the case of 
HAT-P-12b. The third model is a simple flat line that represents 
an atmosphere dominated by clouds that cause gray absorption. 
We multiply all models by the filter transmission curves and 
a theoretical spectrum of HAT -P-12 computed with th e SPEC¬ 
TRUM spectral synthesis code dGrav & Corballvl 19941) to obtain 
theoretical values for k. The match between these three models 
and our data was computed using a chf test. All three models 
were offset in the vertical direction until their mean from the B 
to the z’ band equaled the global k value of Table|3j k — 0.13779, 
see Figure |6] In this way, we set the pressure level that defines 
the planetary radius. The achieved values areXied = 1-61 for the 
flat line model, 2.20 for the cloud-free model, and 5.92 for the 
Rayleigh slope with the number of degrees of freedom equaling 
nine. The corresponding one-tailed probabilities P of a test 
are 0.11, 0.019, and < 1 x 10“^. This indicates that the mea¬ 
sured data points agree with the flat model spectrum, marginally 
agree with the cloud-free model spectrum, but disagree with the 
Rayleigh slope. Moreover, the discrepancy of the g’ and z’ band 
measurements of 3.5 cr points toward a wavelength dependence 
of the planetary radius. 

However, individual data points might still be affected by 
correlated noise. Therefore, we base our results on the ensem¬ 
ble of data points in which individual systematic effects are 
minimized. As a test, we rejected the most deviant data point 
(that of the z’ band) and compared this new spectrum with the 
same three models. After a newly calculated vertical shift, the 
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Fig. 6. Broad-band transmission spectrum of HAT-P-12b. The measured 
values are given in black. Over plotted are a c l oud-fr ee solar metallicity 
spectrum of HAT-P-12b from iFortnev et ak I d201(]h (red solid line), a 
Rayleigh-scattering slope rif?p/riln/l = -4//(blue solid line), and a flat 
line presenting wavelength-independent atmospheric absorption (green 
solid line). When folded with the filter, transmission curves, and a the¬ 
oretical spectrum of HAT-P-12, the model values are given in their cor¬ 
responding color as open circles. At the bottom the transmission curves 
of the used filters are shown. 


^red for rhs cloud-free, and Rayleigh models are now 
0.50, 0.62, and 2.73, respectively. Thus, the data sample is in 
very good agreement with the flat line, but cannot clearly distin¬ 
guish between the constant planetary radius and the cloud-free 
model. However, the corresponding probability for the Rayleigh 
slope of 5 X 10“^ still rules out this model after discarding the 
z’ band. This probability rises slightly if we take into account 
that a constant Rayleigh slope is unlikely to be sustained ove r 
many scale heights owing to grain settling (iPont et al.ll2013h . 
In the case of HD 18973 3b, the Rayleigh slope is constant from 
the UV to the z’ band (ISing et al.ll201 ih . but weakens at red¬ 
der wavelengt hs until k reach es a constant value at about the 
H or K band (iPont et al.ll2013h . To simulate this effect, we ad¬ 
justed the same value of k for the J band as for the F972N20 
filter for the Rayleigh model, which changesto 2.45 and in¬ 
creases the probability to 1 x 10“^. Fixing u to its theoretical value 
reduces the agreement of the measurements with the Rayleigh 
slope model since all data points remain nearly unchanged ex¬ 
cept the point at J, which increases by about 1 cr. A Rayleigh 
slope confined to a narrower wavelength range in the blue fol¬ 
lowed either by a flat line or a cloud-free spectrum at redder 
wavelengths cannot be ruled out by our data. 

In a previous study, iLine et aP (1201 3h obtained a transmis¬ 
sion spectrum from 1.1 to 1.7 pm, which clearly rules out the 
presence of the broad water absorption feature at 1.4 pm pre¬ 
dicted by a cloud-free atmosphere model. They conclude that 
clouds in the atmosphere are the physically most plausible ex¬ 
planation for their rather flat spectrum. Our work extends the 
observed wavelength coverage of HAT-P-12b toward the near 
UV and leaves the cloudy atmosphere as the best explanation for 
the entire spectrum. 

HAT-P-12b has an equilibrium temperature of about 1000 
K. Theoretical models of hot-Jupiter atmospheres predict cloud- 
free atmospheres as clouds of, for example, MgSiOs are ex¬ 
pected to form at hig her pressure levels than th ose probed by 
these observations (see lMadhusudhan et al.l20f4l and references 
therein). However, HAT-P-12b is not the only hot Jupiter of this 


temperature regime that does not fulfill the theor etical predic¬ 
tions. The works on tr ansmission spectrosc opy of iGibson et al.l 
JMl on WAp-29 b, iNikolov et al.l JMl on HAT-P-lb, and 
iMallonn et alJ (1201 5h on HAT-P-19b disfavor clear atmospheres. 
Although for each case there might be individual explanations 
including for example alkali metal depletion by photoionization 
or night-side condensation, one scenario eases the interpretation 
of all these observations: the slant viewing geometry of this ob¬ 
serving technique through the planetary atmosphere strongly en¬ 
hances the opacities of minor condensates. Even in low concen- 
trations, such material can mask atmospheric spectral features 
(lFortnevl2005h . Therefore, in the context of current instrumental 
limits, transmission spectra with muted features might outweigh 
the detections of broad cloud-free features. 

6.2. Periodicity in the monitoring iight curves 

The V band light curve of season 2014 presented in Section 
|5] shows a periodicity of ~ 60 days, which the I band one of 
~ 30 days, thus the question arises if these are linked to the stellar 
rotation period. The semi-amplitudes of the order 1-2 mmag are 
well within the regime of systematic effects. To verify this signal, 
we investigated independent photometric data sets. We down¬ 
loaded the monitoring data of HAT-P-12 of the Super WASP sur¬ 
vey from the NASA Exoplanet Archival observed in 2007 and 
found a periodicity of ~ 65 days, here with a slightly larger am¬ 
plitude of about 3 mmag. However, not far from this period, 
the window function of this data set has a moderate peak at 
~ 80 days. We also analyzed the discovery light curve from the 
HATnet survey, taken in 2006 and another HATnet monitoring 
of this object from the year 2010, both made available on the 
HATnet homepag^l Their photometric precision only rules out 
significant periods with amplitudes higher than about 4 mmag, 
but does not allow for any statement about the potential periods 
at 30 and 60 days. Eurther observations are needed to confirm 
that these potential periods can be attributed to the stellar rota¬ 
tion. 

7. Conclusion 

The scientific goal of this work has been to investigate the trans¬ 
mission spectrum of HAT-P-12b for spectral variations from the 
near-UV (sloan u’) to the near-lR (Johnson J) using broad-band 
filters. We observed 20 transit events with seven different facil¬ 
ities. Two events were observed in multichannel mode, thus in 
total we acquired 23 new photometric transit light curves. We 
also re-analyzed eight transit light curves already published in 
the literature to complete our filter sample, which finally con¬ 
sisted of ten different filters. In addition, the host star HAT-P-12 
was monitored in the observing seasons 2012 and 2014 to eval¬ 
uate its activity level and the level of starspot influence on the 
derived transit parameters. 

The entire sample of transit light curves was homogeneously 
analyzed with the transit modeling software JKTEBOP. The de¬ 
rived transit parameters are generally in good agreement with 
previously published values, but with increased accuracy. Eur- 
thermore, the large number of light curves and their time sep¬ 
aration (2007 - 2014) allows for a refinement of the orbital 
ephemeris. No hint of a periodic transit timing variation was 
found. However, the scatter of the data points in the O-C di¬ 
agram is slightly greater than their mean error bar. Eurther in- 

^ exoplanetarchive.ipac.caltech.edu/docs/SuperWASPMission.html 
^ hatnet.org/planets/ discovery- hatlcs.html 
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vestigations are needed to verify whether the scatter is increased 
by underestimated correlated noise in the data or by an astro- 
physical effect. 

The monitoring data show in general a photometrically quiet 
object. We measured a maximum variability of one percent with 
no periodicity. The data allowed us to exclude any periodic vari¬ 
ations of an amplitude greater than 2 mmag for the monitored 
time span. We concluded that a potential influence of starspots 
on the transit parameters is within the derived transit parameter 
uncertainties, therefore we did not apply an additional starspot 
correction. Our monitoring campaign, together with the HATnet 
and Super WASP photometry of HAT-P-12, yields indications of 
a photometric variation of 60-day period with less than 2 mmag 
amplitude, however the signal is not significant enough to be at¬ 
tributed to the stellar r otation. 

In agreement with iLendl et ^ (1201 3h . we found that multi¬ 
ple transit observations lower the effect of uncorrected correlated 
noise on the derived transit parameters. Because this noise com¬ 
ponent is not correlated between individual light curves, its effect 
is minimized when averaged. Furthermore, the error estimation 
results in more realistic error bars if multiple light curves are an¬ 
alyzed simultaneously, indicated here by the better agreement of 
the transit parameters derived after combining the light curves 
per filter compared to those extracted from the individual light 
curves. 

The combined ten data points over wavelength, which form 
our transmission spectrum, were compared to a flat line, a 
Rayleigh scattering slope (power law coefficient a = 4), and 
a cloud-free solar-composition model. In general, the data favor 
the flat model. However, individual data points deviate signifi¬ 
cantly from each other, especially the z’ band, which shows a 
deeper transit. While additional observations need to confirm the 
reality of this measurement, a Rayleigh scattering model extend¬ 
ing into the NIR can be robustly excluded independently of this 
individual filter. However, a Rayleigh slope confined to shorter 
wavelengths c annot be ruled ou t. Our spectrum strengthens the 
conclusion of iLine et al.l (1201 3h obtained in the NIR that gray 
absorption caused by clouds is the best explanation for the at¬ 
mosphere of HAT-P-12b, and we extend this result to a broader 
wavelength range from the near-UV to the NIR. 
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